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Abstract 

We introduce the ideal Gaussian glass-forming system as a model to describe the thermodynamics 
and dynamics of supercooled liquids on a local scale in terms of the properties of the potential 
energy landscape (PEL). The first ingredient is the Gaussian distribution of inherent structures, 
the second a specific relation between energy and mobility. This model is compatible with general 
considerations as well as with several computer simulations on atomic computer glass- formers. 
Important observables such as diffusion constants, structural relaxation times and kinetic as well as 
thermodynamic fragilities can be calculated analytically. In this way it becomes possible to identify 
a relevant PEL parameter determining the kinetic fragility. Several experimental observations can 
be reproduced. The remaining discrepancies to the experiment can be qualitatively traced back to 
the difference between small and large systems. 

PACS numbers: 64.70.Pf 
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I. INTRODUCTION 



The^iinderstanding of the dynamics of supercooled liquids is still far from being complete 

e, in real space 



mm 

. A lot of insight has been gained from simulations . For exam 
the microscopic nature of dynamic heterogeneities has been clarified 
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Using the framework of the potentia' 



gained in configuration space 
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energy landscape (PEL) a lot of insight could be also 
13| . A key aspect is the use of inherent structures (IS) , 
i.e. local minima of the PEL. Upon minimization basically all configurations can be mapped 
on a IS. In this way the regular dynamics can be mapped on a hopping dynamics between IS 
[l^. [isl. Physically, this mapping can be interpreted as a removal of the vibrational degrees 



of freedom. However, as explicitly shown in 



161 ] the properties of the structural relaxation 



remain identical for sufficiently low temperatures. Generally speaking, the mapping on the 
IS can be interpreted as a coarse-graining procedure. At low temperatures the IS dynamics 
displays many correlated forward-backward jumps between adjacent IS. In a further coarse- 
graining step it is p ossible to define metabasins (MB) by an appropriate merging of adjacent 



IS 



19 



2l\ . In this way the effect of correlated forward-backward motion has 



basically disappeared. 

A key question deals with the relation between thermodynamics and dynamics. For 



example the empirical Adam-Gibbs relation 



22| 



T{T) = Toexp{-BAG/Ts,{T)) 



(1) 



relates the configurational entropy Sc to the local relaxation rate P. A further relation 



between thermodynamics and dynamics is 



brmulated via the fragilities. In the spirit of the 



thermodynamic fragility as discussed in 123 
index via 25 ] 

'^thermo 13 g 



24| one can define the thermodynamic fragility 



(2) 



where Tg = l/j3g (choosing ks = 1) denotes the glass-transition temperature. Furthermore, 
the kinetic fragility is defined as 



rrikv 



dlnTjd{Tg/T). 



(3) 



Qua' 



itatively, it denotes the slope of the relaxation time (or viscosity) in the Angell-plot 



26| . Empirically, one finds a significant correlation between the kinetic and the ther- 



modynamic fragility 



23|. In principle the kinetic fragility may also be defined for the 



27| minor varia- 



diffusion constant. Due to the violation of the Stokes- Einstein relation 
tions of the value of rrikin will be present. Furthermore, it turns out that for the set of all 
glass-forming systems one observes a significant correlation between mfcj„ and the degree 
of non-exponentiality, expressed, e.g., by the exponent Pkww of the stretched exponential 
function 28|]. If one restricts oneself, however, to the set of all molecular glass-forming 



systems (excluding in particular network forming systems and polymers) the residual cor- 
relation is very weak (-0.28) and the values of Pkww are restricted for most of the systems 
(> 80%) in that work to a relatively small regime between 0.5 and 0.62 28|, |29||. In contrast, 
the network-forming systems are characterized by nearly exponential relaxation and small 
values of rrikin. 

In the language of the IS or the MB the thermodynamic properties at constant volume are 



to large extent determined by their energy distribution G{e). For many systems it 



las been 



shown numerically that the distribution of IS can be described by a Gaussian [30|, l31|, l32|, |33 |. 
Even for BKS-Si02 the distribution is Gaussian, albeit displaying a low-energy cutoff in the 
range of accessible temperatures for computer simulations 3J]. Furthermore, it turns out 
that the distribution of IS and MB is nearly identical in the relevant regime of low-energy 



states 



20|. 



Within the PE: 
aspects 



191, 



20 
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approach it is possible to relate the thermodynamic and the dynamic 



35| . This is based on the observation that the escape rate from a MB 



can be expressed in terms of its energy, i.e. r(e,T). Furthermore, it turns out that the 
temperature dependence of the diffusion constant D{T) can be fully expressed in terms of 
the average local escape rate. As a consequence, knowledge of G{e) and r(e,T) allows one 
to predict D{T). A similar type of relation between energy and mobility can be found, e.g., 
for the trap model [sol]. 

The goal of this work is to elucidate the properties of a system with a Gaussian dis- 
tribution G{e) of MB. The functional form of r(e,T) is rationalized by different models, 
discussed in literature, and at the same time by comparison with previous computer simula- 
tions on the binary mixture Lennard- Jones system (BMLJ) and silica (BKS-Si02). On this 
basis we define an ideal Gaussian glass-former (IGGF). For the IGGF several observables 
can be determined analytically such as the temperature dependent diffusion constant and 
relaxation time, its kinetic and thermodynamic fragility and its non-exponentiality. In this 
way it becomes possible, e.g., to elucidate the relevant PEL parameters which determine the 
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fragility. In Sect. II the IGGS is introduced and in Sect. Ill its main properties are calculated. 
We end with a critical discussion and a summary in Sect. IV. 



II. DESCRIPTION OF THE IDEAL GAUSSIAN GLASS-FORMER 
A. Thermodynamics 

Of crucial importance for the properties of a glass-forming system is the number density 
of IS, denoted G{e). Here we always consider a system with particles. For many different 



systems, studied via computer simulations, a Gaussian density of IS has been found [32|, |33|, 



37 



38| 



I.e. 



G(e) = exp(aA^) 



V2 



exp(-(e-eo)V2fx2). 



(4) 



A notable exception is BKS-Si02. This system is characterized by a low-energy cutoff 39 1 
which gives rise to the fragile-to-strong crossover In principle, for the calculations, 

shown below, the effect of a low-energy cutoff can be incorporated 29|]. Here we mainly 
concentrate just on the case of a purely Gaussian density of IS. 

For a closer discussion one has to take into account that the average curvature around the 
minima may depend on e. For different systems it turns out to a very good approximation 
that one has a linear energy-de pen dence for the free energy -Fharm(e), related to the harmonic 



vibration in a well 
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32 
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4J]. This can be written as 



Fharmifi) COTlSt Pharm^^- 



(5) 



The constant Pharm is a material constant. The meaning of the sign of Pharm is visualized 
in FiglH 

The Boltzmann distribution Peq{e) describes the probability to be (at a randomly chosen 
time) in an IS with energy e. Peq{e) is proportional to G{e) exp(— /3e) when Pharm = 0. 
Taking into account the curvature-effect, introduced above, one finds 



Pegie) oc Geff{e) exp(-/3e) 



(6) 



with the effective density 



Geff{e) oc G{e) exp(-Ffe„^^(e)) oc exp(-(e - eo,e//)V2o'^^ 



(7) 
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Pharm < Pharm > 







Figure 1: A sketch of the cases Pharm < and Pharm > 0. Shown are typical curvatures around 
representative IS at different energies. 

and eo,e// = cq + PharmO''^- Thus the presence of an energy-dependent average curvature can 
be incorporated by a shift of the Gaussian distribution of states. 

The standard definition of the configurational entropy is — ^jPilnpj where the sum is 
over all states (not energies). Mapping this relation to the description in terms of energies 
one obtains 

Sc(T) = J depeg{e)\n{G{e)/peg{e)) = j depeq{e)Sc{e) - j depeqie)\npegie). (8) 

For the G{e) and Peq{e), obtained for the Gaussian distribution, one obtains from the first 
term 

S,{T) = Na- {l/2)a^ {(3 - Pharmf ■ (9) 

For large one expects that a'^ (x N due to the central limit theorem. Then Sc{T) becomes 
extensive as expected. In contrast, the last term in EqJHl which would give rise to 1/2, 
can be neglected because it is not extensive and just gives rise to a minor redefinition of a 
(a ^a + l/2N). 

Defining the Kauzmann-temperature by the condition Sc(Tk) = (and Pk = ^/Tr) Eqj9] 
can be equivalently expressed as 



T5,(T) = [(Na) + a'PPK/2 - a'Pi^Um - Tk) 



(10) 



Neglecting the temperature-dependence of the second term this is actually the standard 
expression when deriving the VFT-temperature dependence (e.g. hi{D{T) I Dn) oc 1/(T — Tq) 
where often Tq ~ is found) from the Adam-Gibbs expression 4l(|. Using a similar way 
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Figure 2: Sketch of the multi-step escape process, including the definition of Vq. The barrier with 
the star is supposed to be the critical barrier beyond which pback < 0.5. 



of rewriting the configurational entropy, this type of argument can be already found in 4l| 



In any event, for the further analysis we will use the expression Eql9]due to its simplicity. 
B. Transitions between MB: Models 

There is a long history of models which describe the dynamics in configuration space on 
a phenomenological level 



36 



45 



46 



47 



481, 



49| . One considers a region of the viscous fiuid 



which can cooperatively rearrange via a transition state. For the time being the initial and 
final states may be characterized by the energy of the respective IS (or MB). For sufficiently 
low temperatures the elementary rearrangement process is considered to be activated: the 
system leaves a state with energy e, crosses a high-energy transition state with rate r(e) 
(from now on the variable T is omitted) and ends up in a new state which is uncorrelated 
to the initial one. Different names can be found for essentially identical models (e.g. trap 
model, free energy model) following this scenario. 

The hopping rate r(e) is characterized by two energies. Ccross denotes the energy of the 
IS just after the final barrier, which has a height Vq; see Figj2] for the sketch. According 
to the model assumptions Ccross and Vq are independent of the initial energy e. Actually, 
even in more complex systems like the random energy model one can argue via percolation 
arguments that Ccross is independent of e fsol. More generally, in a percolation-like picture 
of the PEL Ccross corresponds to the energy level from which on the system finds adjacent 
states with similar energies and thus does not have to increase further in the PEL for the 
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final relaxation. Defining Eapp{e) as the apparent activation energy to escape from energy 
e, this scenario can be written as 



T{e) = To{e)exp{-PEapp{e)) (11) 

with 

Eapp{e) = Ccross + Vq - 6 (12) 

for e < Ccross and Eapp{e) = Vq for e > Ccross- Stated differently, the escape for energies lower 
than Ccross is solid-like (activated) whereas otherwise it is liquid-like 35|. 

The energy-dependent prefactor ro(e) reflects possible entropic effects. As argued in 



45 



511] the prefactor ro(e) contains an energy- dependent factor Mentro which denotes the 
number of escape paths to reach a high-energy state with energy Ccross- In most models 
this is neglected by simply choosing ro(e) = Tq. This would be justified in case of ID 
reaction paths or low-dimensional percolation paths. A simple expression for Mentro can 
be formulated if every state with energy Ccross can be reached from exactly one state with 
energy e(< Ccross)- It is given by Mentro = G{ecross)/G{e) Q, i.e. 

ro(e) = ToG{ecross)/G{e). (13) 

This holds for e < Ccross in the opposite limit one just has ro(e) = Tq. For Ccross — e -C 
Co — €.cross Eqil3] can be approximated as 

ro(e) ^ roexp((eo - ecross){ecross - e)/a^). (14) 
For later purposes this is rewritten as 

ro(e) = Toexp{Kkentroiecross " &)) (15) 

with 

kentro \^0,eff (Across 

)/a' (16) 

and 

Co (Across 



K 



cross (17) 
CO,e// Across 

This somewhat complicated way to rewrite Eq{T3]is motivated in two ways. First, because 
Peq{e) is directly related to Geff{e) and thus to eo,e//, see EqJHl it is more convenient to use 
eo,e// rather than cq. Second, in practice the factor k has to be treated as an empirical 




Figure 3: Sketch of the energies, introduced in the text. A possible difference between G{e) and 
G(.ff{e) is neglected. 

parameter because the increase of the entropic term ro(e) with decreasing energy may 
somewhat deviate from the specific scenario, described above. The relevant energy scales 
are summarized in FigJ31 

It is convenient to introduce the shifted inverse temperature 

(3 = (3- KkentTO- (18) 

In principle all calculations, shown in this work, can be performed as well for k ^ 1. However, 
since the influence of the entropic prefactor is not as important as the energetic term the 
additional complexity of the expressions is not worth the additional information for n ^ 1. 
In what follows we therefore always choose k = \. 
Now one can rewrite Eq{TT] as 

r(e) = roexp(-/3(ecross - e)) exp(-/5Vo) 

and, for a Gaussian density of states, 

Peq{e) oc Geff{e) exp(-/5e) oc exp[-(e - e^ross + cr'^Pf /2cr'^]. (20) 

When comparing Eq. [19] with simulations one has to take into account that the simulated 
system may contain more than one elementary system. Each subsystem is characterized by 
an energy Ci and e = ^Cj. For a superposition of M independent subsystems the total 
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(19) 



hopping rate Fa/ (e) is just the sum of the individual hopping rates r(ej). To a first approx- 
imation one may assume that the energy e is equally distributed among the M subsystems, 

lat apart from another energy- 



291]. This expression for rj\/(e) 



yielding Tm^s) = MT{e/M). A closer analysis shows t 
independent factor this is indeed the correct expression 
suggests to generalize EqJT^to 

r(e) = To exp{-Xf3{ecross - e)) exp{-(3Vo) (21) 

Here 1/A is a measure for the number of elementary subsystems, present in the specific 
system. This completes the definition of the IGGF. For later purposes we introduce the 
dimensionless quantity 

/i = /3cr (22) 
which will turn out to be the central quantity characterizing the properties of the IGGF. 

C. Comparison with simulations 

The above scenario has been compared with simulations of relatively small systems of 



the BMLJ system (A^ = 65) [20( and BKS-SiOa (A^ = 99) (52|). This comparison has been 
performed for MB in order to have a random-walk type dynamics in configuration space. 
With IS it would have been impossible to express observables such as the diffusion constant 



or the relaxation time just in terms of the waiting times [20|, |29|. For the comparison the 
average waiting time in MB of a given energy e have been determined, denoted as (r(e)). 
Naturally, the average escape rate r(e) is then given by 

Note that this definition does not imply that the escape from a MB with energy e corresponds 
to an exponential waiting time distribution with average waiting time (T(e)). 

The simulations have fully confirmed the validity of Eql2T] except for a slight smearing 
out effect for energies close to Ccross- Actually, the effective barriers could be identified by 



a closer analysis of the relevant minima and saddles of the PEL [20|. Actually, in [53] it 
has been shown that the additional barrier before the final transition (denoted Vq above) 
and the barrier, governing the local forward-backward motion at low temperatures (within 
a MB) are roughly the same. 







thermodynamic 


dynamic 




N 


a 


^cut 


Pharm 


a 


(Across 


A 




Vo 


To 


BKS-SiOa 


99 


3.5 eV 


43.4 eV 


w 


1.14 


37.5 eV 


0.66 


0.62 


0.8 eV 


1/(20 fs) 


BMLJ 


65 


3.0 




-0.3 


0.73 


12.9 


0.55 


0.3 


1.0 


1/150 



Table I: The thermodynamic and dynamic PEL parameters, obtained from simulation of BKS-Si02 
and BMLJ. 



Very recently, de Souza and Wales have analyzed the temperature dependence of the mean 



square displacement, evaluated for a fixed time r [53|]. Of course, for very large r this analysis 
recovers the standard diffusion coefficient. For ambient r, which for the lowest temperatures 
is significantly shorter than Tq,, the authors observe a simple Arrhenius behavior with the 
high-temperature activation energy Vq. For lower temperatures this approach is sensitive 
to the local forward-backward motion within a MB. The barriers in this regime are of the 
order of Vq so that the local processes remain activated with the high-temperature activation 
energy. This strengthens the observation that it is roughly the same value Vq which governs 
the additional barrier height at low and high energies. 

Furthermore it turns out that ro(e) indeed shows an exponential dependence of energy. 
Interestingly, To^ecross) ~ 1/20 fs~^ is of the order of typical molecular time scales. This 
also suggests that the increase below Ccross is due to entropic reasons. 

The PEL parameters, obtained from the fitting, are listed in TablH Note that if not 
mentioned otherwise from now on all energies are expressed relative to eo,e//; i-e. the maxi- 
mum of Geff{e). For the analytical calculations, to be presented below, it is convenient to 
exclusively use Eqj2Tl i.e. using e < Ccross and P > 0. The first relation starts to be very well 
fulfilled if Ccross ~ {^{T)) > cf which roughly implies T < 0.6 in case of BMLJ and T < 3600 
K in case of BKS-Si02. In this temperature range one also has /3 > 0. 

Interestingly, Ccross is significantly smaller than e^^^ff- As will become clear below this 
difference is crucial for properties like the fragility. The additional barrier height Vq is present 
both for BKS-Si02 and BMLJ (and has similar height after normalization by a). Therefore 
Vo cannot be of any relevance for the question of fragility. It can be directly extracted from 
the high-temperature behavior. 

The observation A < 1 suggests than even these small systems are not elementary. This 
is equivalent to the result reported in 211] that a consistent mapping on an elementary trap 
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model is not possible. 

Two major differences are evident when comparing BKS-Si02 and BMLJ. First, the 
low-energy cutoff for BKS-Si02 is significantly larger than the cutoff, dictated by entropy. 
Thus, the amorphous ground-state is a finite-entropy state. Second, {—ecross) / VAo"^ is much 
lower for BKS-Si02. This means that activated processes become relevant only for states 
much lower in the PEL. As a consequence, a characteristic temperature like Tmct should be 
lower for sihca than for BMLJ. Indeed, ^{g/Tmct) = {cr/TMCT)siiica - {(y /Tmct)bmlj 



12.2 - 6.7 = 5.5 



54 



55| and A{— e cross) / = 6.4 are similar. Furthermore, the energy- 



dependence of ro(e) for BKS-Si02 is much more prominent. 



III. THE DYNAMICS OF IDEAL GAUSSIAN GLASS-FORMING SYSTEMS 
A. General 

n 

The MB dynamics can be characterized by a waiting time distribution ^{t) [20]. From 
this one can calculate the different moments (r") of (fiir). It has been shown in previous 

n 

work that the diffusion constant D is proportional to 1/ (r) [19| . Within the continuous-time 
random walk (CTRW) formalism the structural relaxation time r„ can be identified with 

^ nn 

(r^)/ (r) [56]. Actually, very recently it has been shown |29|,l57|] that it is indeed fully justified 
to use the CTRW-formalism to describe the dynamics of the BMLJ (A^ = 65) system. 

Given the distribution of energies as well as the relation between energy and mobility one 
may ask whether one can explicitly calculate (t"). For this purpose we first introduce (p{e) 
as the probability density that in a series of different MB, visited by the system, a randomly 
chosen MB has energy e. Then the average waiting time is given by averaging (T(e)) over 
all MB, i.e. 

(r) = j deif{e){r{e)) = J de^(e)/r(e). (24) 

(p{e) is distinctly different from the Boltzmann distribution Peq{e) which denotes that at 
a randomly given time the present MB has energy e, i.e. Peqi^) oc V9(e)(r(e)). Including a 
normalization factor this can be rewritten as 

P.(e) - (25) 
Qualitatively, this relation expresses that low-energy states (small r(e)) are often observed 
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(at randomly chosen times) although their actual number oc ip{e) may be very small. Mul- 
tiplication of Eql25]with r(e) and subsequent integration yields 

(r)-i = J depe,{e)Tie) ^ {T),. (26) 

Thus, the average waiting time is also related to the rate average over the equilibrium 
probability distribution. Note the different notations ((.) as the yj-average vs. {.)p as the 
p-average. Using the explicit form of Geff{e) one obtains after a straightforward integration 

(r)-i = roexp((AV2 - A)/iV2) exp{-f3Vo) (27) 

So far no information about the nature of the relaxation process has entered the analysis. 
In the most simple case the escape from a state with energy e is governed by a single barrier 
height. Then the waiting time distribution, related to this energy, is just r(e) exp(— r(e)t). 
For the BMLJ(N=65) system one has 1/A ~ 2 subsystems. In the most simple picture the 
total energy is then the sum of two independent subsystems, each with energy Cj (ei + 62 = e) 
and for a given energy decomposition the total rate r(e) is given by r(ei) + r(e2). Actually, 



as outlined in 29|, the normalized second moment (T(e)^)/(r(e))^ is expected to be around 
16 for T = 0.5 for 2 subsystems as compared to 2 for an elementary system. The broadening 
of the waiting time distribution at fixed energy is due to the fact that for a given total 



energy e several decompositions e = ei + 62 are possib 



e, each giving rise to different escape 



rates. The numerically observed value is approx. 8 [58|]. This means that the BMLJ(N=65) 
system behaves, to first approximation, like two independent subsystems (each described by 
A = 1 and variance cr^/2 if cr^ is the variance of the original system). A possible reason for 
the decrease of 16 to 8 will be given below. In any event, in what follows we neglect this 
effect and postulate that the elementary system behaves like an IGGF with A = 1 and an 
exponential waiting time distribution at given energy. Since the waiting time distribution 
at fixed energy is a well-defined observable in the MB approach the subsequent calculations 
could be easily generalized to take into account deviations from a purely exponential behavior 
of the waiting time distribution of the elementary system. 

This aspect is strongly related with the old discussion of homogeneous vs. heterogeneous 



relaxation 59 



60l |. Heterogeneous relaxation would simply mean that one has a superposi- 



tion of exponentially relaxing entities. Experimental 



at the glass transition is basically heterogeneous 6l|]. This indicates that the choice of an 



y it has been shown that the dynamics 



exponential waiting time distribution is indeed not too bad. 
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B. Calculation of moments 



With this approximation the waiting time distribution y9(r) and the distribution v?(e), 
reflecting the thermodynamics, are related via 

(^(r) oc j de j dt Lp{e) exp{-T{e)t)5{t - r). (28) 

Its different moments (r") can be directly calculated 

(r") = [ deip{e)n\T{e)-'' = n!(r)(r^-")pexp(n/5\/o). (29) 



For the second equality Eql25] has been employed. 

Straightforward evaluation of Gaussian integrals yields 

{{T/ToDp = exp[(mV2 - m)//^ - m(3Vo]. (30) 

The case m = 1 recovers EqlTTl Furthermore, the case m = — 1 gives finally rise to 

(r2)/(r)2 = exp(/.2) (31) 

In most models no distinction between Ccross and eo,e// is made. Then /i can be identified 
with f5. The relations for this special case can be already found in literature . Note that 
in this limit EqJTTl corresponds to the well-known temperature dependence, discussed, 
e.g., in [50 1 . 



IV. APPLICATIONS 
A. Kinetic fragility 

Here we analyse the temperature dependence of (r) (and thus of D{T)) and in particular 
the fragility. The glass transition temperature is defined by the condition 

ro{r{T,,K)) = 10^. (32) 

Neglecting for the time being the somewhat different temperature behavior of D{T) and f]{T) 
(see below) T^^ie = l//5g,i6 roughly corresponds to the calorimetric Tg because rj{Tg)/rj{T ^ 
Tg) ^ 10^^. Simple expressions emerge for the case Vq = (corrections can be simply 
calculated but only mildly influence the results). Using Eql30]a simple calculation yields 

(^l3g,K = kentroo + sj2K\n{lQ). (33) 
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Figure 4: The temperature dependence of ro(T)(oc D{T)) for different values of the crossover 
energy with A = 1 (the values are given with respect to eo,e//)- 

In relation to the definition of Tg we use the notion mkin,K rather than nikin (see Eqj3]) to 
express the dependence on the time scale. Then a straightforward calculation yields 



2K + ^y2K/\n{W)k,ntroCr. 



(34) 



In this regime the fragility depends on the dimensionless parameter kentro<^ = ~(^crossl<^- 
Thus, the dynamic crossover energy is a central PEL parameter determining the fragility. 
These results are visualized in FigjH One can clearly see how the fragility increases with 
increasing — Ccross/o"- 

Note that EqJM] implies that BML J is stronger than BKS-Si02 if the cutoff were artifi- 
cially removed so that the PEL is purely Gaussian. The non-fragile behavior of BMLJ has 



been already mentioned in Ref. 



62| 



Of course, since the temperature dependence of is in general not identical to that of 
(r) the results would slightly differ if mkin,K is calculated for Tq, or r] rather than for the 
diffusivity. 

Empirical relations to correlate the fragility with, e.g., the Poisson ratio have been sug- 



gested 



631 ] but are questioned in 



6J]. It would be interesting to check whether there exists a 
physical connection between the observables, suggested in that work, and the value of Ccross, 
determining the crossover from liquid-like behavior to solid-like behavior. 
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B. Relation to the AG approach 



Alternatively, one can calculate the value of (3g under the assumption of the AG relation 
EqUJand a Gaussian PEL (using jSharm = 0). Then one has to solve the equation 



10^ = exp{PgBAG/{a - Py/2N). (35) 



For large K one obtains 



V2aN BagVN 



Then a straightforward calculation yields for the fragility (again in the limit of large K) 



_ v^ir^(lnlO)V 

mkin,K — 7= • [<ii) 

BagV N 



Within the AG-approach the fragility depends on the density of states, i.e. a, as well as the 
empirical constant Bag- A large number of states implies larger fragility (at least for fixed 
Bag which, of course, could also depend on a [25I). 

It may be interesting to compare this relation with the fragility EqJMj obtained for an 
IGGF. Qualitatively, both relations would show a somewhat similar behavior if systems with 
large a are related to systems with a low crossover energy Ccross, i-e. large kentro- This is not 
unreasonable because in the spirit of percolation-like arguments for a larger number of IS the 
system would be able to find a path with a lower barrier to move between two low-energy IS. 
However, in a strict way it will not be possible to map Eql37]on Eq|34]because of the different 
i^- dependence. Formally, this problem could be solved if a decreases with increasing K, 
i.e. going to longer time scales and thus lower glass transition temperatures. Qualitatively 
this statement is equivalent to the requirement that G{e) decays faster than a Gaussian. 



This has been suggested in [651]. Physically this might, e.g., occur as a consequence of a 
broadened low-energy cutoff. 



C. Thermodynamic fragility 



In the spirit of the thermodynamic fragility as discussed in 23|, |2J] one can define the 



thermodynamic fragility index via 



25| 



'^thermo^K 



-A 



(38) 
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We obtain, using Eqj9], 

Note that the denominator must be positive because otherwise the entropy of the system 
would be negative. Under this condition, an increase of crkentro (which is the only relevant 
dimensionless parameter, characterizing IGGF) and thus of aPg (via Eq. [33]), gives rise to an 
increase of mthermo and rrikin, independent of the values of Pg or a. This strong correlation of 
rrikin and mthermo is in agreement with the experimental observation for most systems 23 1. 

Interestingly, increasing the value of Pharm yields a decrease in mthermo,K- However, a 
different behavior emerges if one includes the vibronic contribution into the entropy, i.e. 
by using SexiT) = SdT) + SharmiT) rather SdT). A straightforward calculations yields 
SexiT, Pharm) = SJT, —Pharm), thereby neglecting a constant and a term, depending loga- 
rithmically on P 291]. Accordingly, when defining rrikin^K on the basis of Sex{T) one obtains 
an increasing thermodynamic fragility for increasing Pharm in agreement with the qualitative 



discussion in 



23|. 



If the cutoff starts to influence the system a detailed calculation is no longer possible 
because the behavior of the configurational entropy at low temperatures depends on the 
details of G{e) at low energies. Thus, it is not surprising that for Si02 the thermodynamic 
fragility does not follow the general trend j23|. 

The present discussion complements the work by 4l(] where the kinetic and the thermo- 
dynamic fragility have been discussed with reference to the AG-relation. Simulations have 
also revealed a significant correlation between both fragilities. 



D. Relaxation properties 



Here we ask for the probability So{t) that a system in equilibrium has not performed a 
hopping process until time t. It is given by 



So{t) = J depeg{e)exp{-T{e)t). (40) 

In what follows the trivial factor exp(/9Vo) will be omitted. For sufficiently low temperatures 
the decay of this function can be related to the structural relaxation 56|, l57 |. 



As shown in 



29 



Q^ one can approximate for intermediate times {So{t) ~ 1/e) 
5o(t) ^exp(-(t/rwH/)^™) 



(41) 
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with 

Pkww = l/VlTJ^ (42) 

and tkww = 1/r* where 

r = r((e(T)))=exp(V)- (43) 

This may justify the use of the stretched exponential as a fitting function at least for 
intermediate times. This result is insensitive to the specific form of r(e) since r(e) only 
enters via T*. Note that for the IGGF the non-exponentiality tends to increase when going 



to 



29 



ower temperatures. Furthermore one can show that in very long-time decay is algebraic 

5o(t) ocr"/''^'. (44) 
One can define the a-relaxation time via 

r„= /" dtSoit) (45) 



which corresponds to the typical time until a particle jumps for the first time 56|, 1571]. From 
Eq. Ho] one immediately obtains (also using Eqj30l) 



Ta = (r-^)p = (l/Fo) exp(3/iV2). (46) 

This has to be compared with the average hopping time (r) (Eq j271) . One obtains 

rj(r)=exp(/i^). (47) 

Since the left side is proportional to Dt^ Eqj47] expresses the invalidation of the Stokes- 
Einstein relation for IGGF. Using the definition of the exponent a via D{T)Ta{T) oc , 
i.e. (r)(r^^) oc (r~^)" one obtains a = 2/3. Experimental values are smaller (e.g. 0.25 for 
orthoterphenyl j27] and 0.23 for TNB 67|). Thus, the decoupling seems to be too strong. 



Qualitatively the strong increase of with decreasing temperature is due to the very long- 
time tail of So(t). 



V. DISCUSSION AND SUMMARY 



The IGGF has been introduced, based on the numerical results for BMLJ and BKS- 
Si02 (except for the low-energy cutoff for BKS-Si02) at small system sizes. The general 
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concepts are also compatible with several models proposed to rationalize the dynamics of 
supercooled liquids. Thus, one naturally finds how properties such as the non-exponentiality 
are generated. 

More specifically, the key conclusions are as follows: 1.) If the cutoff-energy does not 
interfere the temperature-dependence of the dynamics is fully captured by the value of n 
(except for a trivial exp(— /5Vo)-term). This means in particular that at Tg an IGGF has a 
fixed value of /i, independent of crkentro and thus independent of its fragility. This implies via 
Eqj42] that the stretching parameter I3kww does not depend on the fragility if determined 
exactly at Tg. This may rationalize the weak correlation between (3kww and rrikin for the 
molecular glass-forming systems, as mentioned above. Of course, residual fiuctuations are 
expected when the smaller-order effects of A, k and Vq are taken into account. 2.) The 
fragility of a system is to a large extent dominated by the crossover energy Ccross relative to 
the width of the energy distribution, i.e. a. Systems are more fragile if the crossover from 
solid-like activated dynamics to liquid-like non-activated dynamics occurs at low energies, 
relative to the width of G{e). Of course, as soon as the low-energy cutoff of the PEL comes 
into play (such as for BKS-Si02) the system automatically behaves Arrhenius-like and thus 
is classified as strong. This also shows that the fragility is only partly able to classify a glass- 
forming system because already the present discussion shows that there at least two very 
different parameters, Ccut and Ccross^ which strongly influence the fragility. 3.)Although the 
BMLJ data can be flt to the AG-relation, from a conceptual point of view the IGGF is not 
compatible with the AG-relaxation. This can be seen from the different dependence of the 
fragility on K. On a qualitative level this discrepancy could be reduced if the distribution 
of states decays stronger than a Gaussian at the low-energy end. 4.) The thermodynamic 
fragility indeed is correlated with the kinetic fragility, albeit in a non-trivial way. Again, the 
systems with a cutoff-behavior (most notably BKS-Si02) have to be discussed separately. 
5.) Finally, the IGGF displays non-exponential relaxation with a non-exponentiality which 
increases with decreasing temperature and, in agreement with the experiment, shows a 
violation of the Stokes-Einstein relation. 

Conceptually, the presence of individual relaxation processes naturally is attributed to 
small systems, reflecting the typical length scales of cooperative dynamics during single MB 
transitions. Thus, any strict comparison with simulations in the framework of the PEL 
approach is conveniently performed with small systems. As shown in previous work the 
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diffusion constant as well as the thermodynamic properties of the BMLJ (A^ = 65)-system 
only have very minor finite-size effects when comparing with the results obtained for much 



larger systems 3l|, |68|, |69|] . However, the structural relaxation time as the well as the non 



exponentiality has somewhat larger finite-size effects 



69| . This effect can be understood 



if one assumes a specific type of coupling between adjacent subsystems of a larger system. 
When some subsystem relaxes it may charige the mobility of the adjacent subsystems 29 1. 
A similar idea can be already found in [sa, 0] and has been also implemented in the 
context of the rate memory to explain the results of multidimensional NMR experiments 
0,0]- ^his very immobile regions typically become mobile at some stage 



48 
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and can relax subsequently. In some sense this idea is also related to the philosophy of the 
facilitated spin models where the local mobility is also influenced by the state of the neighbor 
spms The coupling between adjacent subsystems can be formulated such that 

the diffusion constant and the thermodynamics does not change whereas the structural 
relaxation time, all moments (r") for n > 2 and the degree of non-exponentiality decrease 
upon this coupling [29]. This might also explain why the second moment for the BMLJ 
system is by a factor of 2 smaller than expected (see above). In particular the exponent 

a, characterizing the violation of the Stokes-Einstein equation approaches experimentally 

i 

relevant values |29i]. However, one of the key results, namely the utmost relevance of a 
single dimensionless parameter fi would still be valid. In any event, the path back from 
small systems to macroscopic systems is one of the challenges for future work. Using the 
IGGS as the elementary system for such models is definitely a reasonable starting point. 

We gratefully acknowledge important input from C. Rehwald and O. Rubner as well as 
very helpful correspondence with L. Berthier about this topic. 
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